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^ ABSTRACT 



The dilepton production process at hadron colliders in the Randall-Sundrum (RS) 
model is studied at next- to-leading order in QCD. The NLO-QCD corrections have 
been computed for the virtual graviton exchange process in the RS model, in addi- 
tion to the usual 7, Z-mediated processes of standard Drell-Yan. K-f actors for the 
cross-sections at the LHC and Tevatron for differential in the invariant mass, Q, and 
the rapidity, Y, of the lepton pair are presented. We find the K-f actors are large over 
substantial regions of the phase space. 
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In brane-world models, the four dimensional universe is a dynamical hypersurface: 
a D 3 -brane (or 3-brane) existing in a higher dimensional spacetime. In many such 
models, the Standard Model (SM) fields are localized on the brane and only gravity 
can propagate in the bulk. The scale of quantum gravity can be lowered down from 
the Planck scale to the TeV scale in these models [1] making it exciting for high- 
energy physics not only because these suggest fresh perspectives to the solution 
of the hierarchy problem but also because these models throw open the possibil- 
ity of the discovery of new physics at energies accessible to collider experiments. 
In addition, these models provide new frameworks for gauge symmetry and super- 
symmetry breaking and suggest theoretical approaches to the cosmological constant 
problem and dark-matter problem. 

The simplest model seeking to address the gauge hierarchy problem was the 
the ADD model proposed by Arkani-Hamed, Dimopoulos and Dvali [2], where, 
starting from a higher dimensional theory, an effective four-dimensional theory at a 
scale M s ~ TeV is obtained. This is done by compactifying the extra dimensions to 
magnitudes which are large compared to the Planck length [3]. 

The main problem that one faces within the ADD model is the reappearance 
of disparate scales viz., the string scale and the inverse of the compactification ra- 
dius. It was an attempt to avoid this problem that led to the formulation of the 
Randall-Sundrum (RS) model [4]. In the RS model the single extra dimension <p is 
compactified on a S 1 /Z 2 orbifold with a radius R c which is somewhat larger than 
the Planck length. Two 3-branes, the Planck brane and the TeV brane, are located 
at the orbifold fixed points = 0, it, with the SM fields localised on the TeV brane. 
The five-dimensional metric, which is non-factorisable or warped is of the form 

ds 2 = e-^r/^dx^dx" + R 2 c d<f) 2 . (1) 

The exponential warp factor e~ KRc<l> serves as a conformal factor for fields localised 
on the brane. Thus the huge ratio ~ 10 15 can be generated by the exponent JCR C 
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which needs to be only of O(10) thereby providing a way of avoiding the hierarchy 
problem. There remains the problems of stabilising R c against quantum fluctuations 
but this can be done by introducing an extra scalar field in the bulk [5, 6]. 

The tower of massive Kaluza-Klein (KK) excitations of the graviton, hffi, interact 
with the SM particles by: 

1 pTTlCRc CO 

An* ~ -j^T^(x)h^(x)-—Y l T^(x)h%\x). (2) 

T^ v is the symmetric energy-momentum tensor for the particles on the 3-brane. The 
masses of the hffi are given by 

M n = x^e^, (3) 

where the the zeros of the Bessel function J\(x). The zero-mode couples 

weakly and decouples but the couplings of the massive RS gravitons are enhanced 
by the exponential e wlCRc leading to interactions of electroweak strength. Conse- 
quently except for the overall warp factor in the RS case, the Feynman rules in the 
RS model are the same as those for the ADD case [7, 8]. 
The basic parameters of the RS model are 

c = K/Mp, (4) 

where m is a scale of the dimension of mass and sets the scale for the masses of 
the KK excitations, and c is an effective coupling. The interaction of massive KK 
gravitons with the SM fields can be written as 

dnt °-J2T^(x)h$(x). (5) 

m n 

Since K, is related to the curvature of the fifth dimension we need to restrict it to 
small enough values to avoid effects of strong curvature. On the other hand /C 
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should not be too small compared to M P because that would reintroduce a hierar- 
chy. These considerations suggest 0.01 < c < 0.1. For our analysis we choose to 
work with the RS parameters c and M 1 the first excited mode of the graviton rather 
then m . 

The decoupling of the graviton zero-mode and the existence of a mass gap in the 
spectrum of KK gravitons imply that it is only the resonant production and decay 
of the heavier KK modes or the virtual effects of the KK modes that one can hope to 
detect in collider experiments. The phenomenology of resonant production of the 
KK excitations and the virtual effects have already been studied in processes like 
dilepton production [9], diphoton production [10], ti production at hadron colliders 
[11], r-production at a linear collider [12] and pair production of KK modes in e + e~ 
and hadron hadron colliders [13]. The sensitivity of the CMS experiment to the 
resonant production of RS graviton KK modes has been studied for electron pair 
production [14]. Recently D0 has reported the first direct search for RS graviton 
KK modes using dielectron, dimuon, and diphoton events [15]. 

In an earlier work we had presented NLO-QCD corrections for e + e — > hadrons 
[16] and dilepton pair production at hadron colliders [17] in the ADD model. These 
results for the dilepton pair production case are extended to the RS model, in this 
paper. We note that it is the same virtual graviton exchange process that contributes 
to dilepton production in both the ADD and RS models. The leading order process 
being the same, the QCD corrections are also not model-dependent. However, as 
explained above, the differences between the two models arise because of the dif- 
ference in the summation over the tower of KK gravitons and also in the overall 
factors. Consequently, the relative weight of the subprocess cross-section due to 
graviton exchange vis-a-vis the SM subprocess will be different in the two models. 
This results in different K-f actors in the ADD and RS models and the dependence of 
the i^-factors on the kinematic variables are also different. In this letter, we present 
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the results for dilepton production at the LHC and Tevatron in the RS model. 

The process we are interested in is where two hadrons Pi , P 2 scatter and give 
rise to leptonic final states, say fi + , [T 

Piipi) + P2(P2) - + *T (Z 2 ) + X(P X ) , (6) 

where p\ , p 2 are the momenta of incoming hadrons P\ and P 2 respectively and fj,~ , /i + 
are the outgoing leptons which have the momenta h,l 2 . The final inclusive hadronic 
state is denoted by X and carries the momentum P x . The hadronic cross section can 
be expressed in terms of partonic cross sections convoluted with appropriate parton 
distribution functions as follows 

2S^^(r,Q 2 ) = £ dxj dx 2 f*( Xl )f»(x 2 ) 

ab=q,q,g 

,1 _ da ab 

X 



J dzls-^ (z, Q 2 ) 5{t - zx x x 2 ) . (7) 



dQ 2 

The scaling variables are defined by k\ = xipi, k 2 = x 2 p 2 where ki, k 2 are the mo- 
menta of incoming partons. 

( Pl +p 2 ) 2 = S, (k l + k 2 ) 2 = s, (h + k) 2 = q-q = Q 2 , 

T = — , Z = — , T = X\X 2 Z . (8) 

b s 

m 

The partonic cross section for the process a(ki)+b(k 2 ) — > j(— g) pi) is given 

i 

by 

2s ^ = J- 2 fdPS^ 1 \J^»y-P j (q)-^ l {q)-C»'^ +l -(q). (9) 

In the above equation, the sum over Lorentz indices between matrix element squared 
and the propagators is implicit through a symbol "dot product". The m + 1 body 
phase space is defined as 



x(27T) n 5^(h + k 2 + q + J2 Pi ) : (10) 

i 

where n is the space-time dimension. The propagators are 



1 



PzM = - (Q--Ml- 1 M Z T Z ) S >'-' (12) 



Po(q) = »W 2 )-S„„A„fe) . (13) 



where 



r ( \ - ( Me\ ( _ q » q A , ( _ mA\ I mA 
~ \ 9flp M n 2 )\ 9va M n 2 ) + \ 9pu M n 2 )\ 9up M 2 ) 

-—i{^-M7){ 9 "-m- (14) 

The function V(Q 2 ) in the graviton propagator Eq. (13), results from summing over 
the KK modes, given by 

V(Q2) - g, w^ ^- (15) 

where M n are the masses of the individual resonances and the r„ are the corre- 
sponding widths. The graviton widths are obtained by calculating their decays into 
final states involving SM particles. A is defined as 

OO v»2 2 ' r rf 

Ju „ •X,' 6 «X > 7 

= E 2 2 ^ r? , (16) 

n=l x s x n mo x « 

where x s = Q/m . We have to sum over all the resonances to get the value of X(x s ). 
This is done numerically and for a given value of x s , we retain all resonances which 
contribute with a significance greater than one per mil, and treat the remaining KK 
modes as virtual particles (in which case the sum can be done analytically). 

We now present the distributions in the invariant lepton pair mass, Q, and the 
rapidity of the lepton pair, Y at the LHC (VS = 14 TeV) and Tevatron (VS = 1.96 



TeV). From these distributions the effects of the NLO-QCD corrections can be clearly 
discerned. For the parton density sets we adopt in leading order (LO) MRST 2001 
LO (A = 0.1670 GeV) and in next-to-leading order the MRST 2001 NLO (A = 0.2390 
GeV). For LHC we choose the kinematic ranges 300 GeV < Q < 3000 GeV and 
\Y\ < 2.2 at Q = 1.5 TeV. For Tevatron 300 GeV < Q < 1000 GeV and \Y\ < 0.9 at 
Q = 300 GeV. The renormalisation scale is taken to be same as the factorization scale 
fi F and /j, F si chosen to be fi F = Q. 

The cross-section dojdQ as a function of Q to NLO is presented in Fig. la for 
LHC. For the figure, we have chosen the representative values of the RS model 
parameters: M 1 = 1.5 TeV the first RS resonance mass and the coupling constant 
c = 0.01. The width of the resonance is related to c and hence a smaller c cor- 
responds to a narrow resonance. The subsequent resonance are determined by m 
and x n . To LO the dilepton case has been presented in [9]. To see the effect of the 
NLO effect we study the fT-factor for the Q and Y distribution. 

The fGfactor for the invariant lepton pair mass distribution defined by 

K = [~dQ~\ [ dQ \ ' (17) 
where I = SM, I = SM + GR for both SM and gravity combined and I = GR 

for only gravity. It is possible to define K GR for the invariant lepton pair mass 

distribution, as there is no interference with SM [17]. The results are presented in 

Fig. lb. The parameters chosen are the same as in Fig. la. In order to understand 

the behaviour of K-i actor of the model involving both SM and gravity, it is useful 

to express it as 

R SM + K GR K (0) 



K (SM + GR) {Q) 



(18) 



1 + KM 

where we have introduced a quantity K^ Q \ defined as the ratio of the LO distribu 
tion of gravity to SM, given by 



K (0 \Q) = 



aa LO 



(Q) 



l-lr 



dQ 



dQ 



(19) 



The behaviour of K 1 ^ (Q) is governed by competing couplings constants of SM and 
gravity and the parton fluxes. In the RS case the gravity contribution is significant 
in the resonance region, (see Fig.la). In the off resonance region the K-factor is 
hence purely K SM . In the resonance region where the gravity effect dominates the 
^(sm+gr) f ac t or shifts to the K GR value (see Fig. lb). This behaviour of the fT-factor 
of the RS case is very distinct from the corresponding case we presented in the ADD 
case [17]. To incorporate the NLO effects for an appropriate distribution one needs 
to take into account the behaviour of the i^-factor accordingly. For Mi = 300 GeV 
the K-{ actor is about 1.5 in the resonance region. This is due to the fact that at loq 
Q (Q = 300 GeV) the gluon flux becomes dominant at Tevatron. The behaviour of 
K GR is the same as in the ADD case [17]. 

In Fig. lc, we have plotted the scale variations of the Q distribution for both LO 
and NLO cross sections. We define R 1 for the invariant lepton mass distribution as 



r lo - 



dvlo(Qil J ' = Vo) 


-i 




dQ 




dQ 



R 1 



NLO 



d<T NLO {Q,V = (M)) 



dQ 



dQ 



(20) 



where /i is a fixed scale which is chosen to be /i = 1.5 TeV for LHC. As can be 
seen from the figure, the inclusion of the NLO corrections stabilises the cross-section 
with respect to the scale /i. Here we have chosen /j, = 1.5 TeV, ie. the first resonance 
region. The scale variation is driven by the gravity part as its the dominant contri- 
bution. 

In Fig. 2a the double differential cross section d 2 a/dQdY is displayed for rapid- 
ity region \Y\ < 2.2 for a Q value of 1.5 TeV. To plot this distribution the Q value 
is chosen such that it lies at the first resonance, where the gravity effect dominates. 
Hence the dominant contribution is purely gravity. The RS model parameters re- 
main the same as before. The K factor as a function of Y is plotted for a choice 
of Q in the resonance region where the dominant contribution to K ( - SM+GR ^ factor 



comes from the gravity part Fig. 2b. The R ratio using the Y distribution is plotted 
in Fig. 2c for the central region of rapidity and for a Q value of 1.5 TeV. In this region 
the scale variation is also dictated by the gravity contribution. 

The corresponding analysis for the Tevatron is done for the Q range 300 < Q < 
1000 GeV and for the RS parameter Mi = 300 GeV and the coupling c = 0.01. 
At low Q the gravity effects of the RS model is dominant in the resonance and off 
the resonance region the effect is negligible. As Q increases the effect of gravity 
starts to become comparable to the SM contribution as is seen towards the third 
resonance in Fig. 3a. In Fig. 3b we have plotted the /^-factor for Q distribution 
at the Tevatron. Using Eq. (18) we can understand the behaviour of K SM+GR . As 
expected the behaviour of K GR is same as the ADD case [17]. The double differential 
cross section d 2 a/dQdY for Q = 300 GeV is plotted as a function of rapidity Y. In the 
resonance region the dominant contribution is from the RS. In contrast for ADD [17] 
only at large Q the gravity effects became comparable to the SM at Tevatron. The 
corresponding ^-factor is plotted in Fig. 4b. Scale variation for the Q and double 
differential dQdY is given in Fig. 3c and Fig. 4c respectively. 

In a recent analysis by D0 [15], the LO cross section was scaled by a constant 
K-i actor of 1.34 to account for the NLO effect for the RS case. This does not yield 
a realistic picture as can be seen from Fig. 3b. In the RS case due to the resonant 
production, the K-i actor is very different from the ADD case reported earlier [17]. 

In summary, we have presented the results for the cross-section for dilepton 
production in the Randall-Sundrum model at the LHC and Tevatron. The large 
incident gluon flux at the LHC makes the NLO QCD corrections very important. 
Moreover, when the NLO corrections are taken into account the cross-sections are 
stabilised with respect to scale variations. In order to derive robust bounds on the 
RS model at the LHC using the dilepton production process, the inclusion of the 
NLO QCD corrections in the cross-section is crucial. 
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Figure Caption 



Figure 1. (a) The cross section is plotted as a function of invariant mass Q of the lep- 
ton pair for Mi — 1.5 TeV at LHC. (b) The corresponding fT-factor for Q distribution 
SM, gravity and SM plus gravity, (c) Scale variation of the cross section at LO and 
NLO as defined in Eq. (19) for Q = 1.5 TeV. 

Figure 2. (a) The double differential cross section d 2 a/dQdY is plotted as a function 
of rapidity Y f or Q = 1.5 TeV at LHC. (b) The K-f actor for the distribution in (a) is 
plotted for the rapidity range, (c) The scale variation of the ratio R is plotted as a 
function of /Va*o for Y — 0. 

Figure 3. (a) The cross section is plotted as a function of invariant mass Q of the 
lepton pair for M 1 = 300 GeV and c = 0.01 at the Tevatron. (b) The fT-factor for 
Q distribution for the same RS parameters in (a) is plotted, (c) The variation of the 
cross section with respect to the scale. 

Figure 4. (a) The double differential cross section d 2 a/dQdY is plotted as a function 
of rapidity Y f or Q = 300 GeV at the Tevatron for the RS parameters M x = 300 GeV 
and c = 0.01. (b) The K-f actor for the distribution in (a) is plotted for the rapidity 
range, (c) The scale variation of the ratio R is plotted as a function of pi/no for the 
central rapidity region Y — 0. 



10 



References 

[1] P. Horava and E. Witten, Nucl. Phys. B460 (1996) 506; J.D. Lykken Phys. Rev. 
D54 (1996) 3697; E. Witten, Nucl. Phys. B471 (1996) 135. 

[2] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B249 (1998) 263; 
I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B436 
(1998) 257. 

[3] For a review of the ADD model, see I. Antoniadis and K. Benakli, Int. J. 
Mod. Phys A15 (2000) 4237; A. Perez-Lorenzana, Univ. of Maryland Preprint 
No. UMD-PP-00-088 (2000) hep-ph/ 0008333. K. Sridhar, Int. }. Mod. Phys A15 
(2000) 2397. 

[4] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370. 

[5] W.D. Goldberger and M.B. Wise, Phys. Rev. Lett. 83 (1999) 4922; Phys.Lett. 
B475 (2000) 275. 

[6] C. Csaki, M. Graesser, L. Randall and J. Terning, Phys. Rev. D 62 (1999) 045015, 
C. Csaki, M. Graesser and G.D. Kribs, Phys. Rev. D 63 (2001) 065002. 

[7] G. E Giudice, R. Rattazzi and J. D. Wells, Nucl. Phys. B544 (1999) 3. 

[8] T. Han, J. D. Lykken and R-J. Zhang, Phys. Rev. D59 (1999) 105006. 

[9] H. Davoudiasl, J.L. Hewett and T.G. Rizzo, Phys. Rev. Lett 84 (2000) 2080; ibid. 
Phys. Rev. D63 (2001) 075004. 

[10] K. Sridhar, JHEP 0105 (2001) 066. 

[11] S. Lola, Prakash Mathews, Sreerup Raychaudhuri and K. Sridhar, CERN 
Preprint No. CERN-TH-2000-275 (2000) hep-ph/ 0010010. 



11 



[12] N. Mahajan, J.Phys.G 29 (2003) 2677. 

[13] Pankaj Jain and Sukanta Panda JHEP 0403 (2004) Oil. 

[14] C. Collard and M.-C. Lemaire, CMS NOTE 2004/024. 

[15] V. M. Abazov et. ah, D0 Collaboration, hep-ex/0505018. 

[16] Prakash Mathews, V. Ravindran and K. Sridhar, JHEP 0408 (2004) 048. 

[17] Prakash Mathews, V. Ravindran, K. Sridhar and W.L. van Neerven Nucl. Phys. 
B713 (2005) 333. 



12 



n 1 r 



10 



-2 - 

10 =~~ 



da/dQ (pb/GeV) (LHC) 



M 1= 1.5 TeV 
c =0.01 




-9 - 

10 E- 



SM 

_ _ GR 
SM+GR 



\ / 



10 



-10 



j i i_ 



1000 



[ 

2000 



L 

3000 



Fig. la 



13 




14 
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Fig. 2b 



17 



1.6 



1 



0.8 



1 1 

R-Ratio (Y) (LHC) 


i i 
i 

M.=1.5 TeV 
i 




c n =0.01 




* 
• 


u n =1.5TeV " 


* 

— • 

• 
* 

— * 


v— n 


• 
• 

— • 




• 




• 

* 

• 




• • « 

— • 


GR LO 


• 


GR NLO 


* • ■ • • 


(SM+GR) LO 


\ \ 

\ \ 
\ m 
\ \ 


(SM+GR) NLO 


\ • 

~ \ \ 
























1 .... 1 . 






0.5 1 1.5 2 

|U/|Ll 
Fig. 2c 



18 



10 

1 

-1 - 
10 r 



10 r 



T 



n 1 r 



T 



~i 1 r 



da/dQ (pb/GeV) (Tevatron) 



M^OO GeV 
c =0.01 




10 



10 



10 



10 r 



10 r 



10 



10 



-9 - 



_I I L 



J 1 l_ 



200 



400 



600 Q 800 



1000 



Fig. 3a 



19 




Fig. 3b 
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Fig. 3c 
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